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ABSTRACT 

The  theory  based  on  covalent  electron  transfer  (CET)  between  the  dx 2_y2  orbitals  of 
mixed-valence  Cu  ions  that  was  introduced  previously  to  explain  the  variation  of  critical 
temperature  Tc  with  chemical  composition  and  the  temperature  dependence  of  high-7^  super¬ 
conductor  properties  is  modified  to  reflect  lattice- vibration  scattering  of  the  large  polarons  that 
dominate  the  electrical  conductivity  of  these  compounds  at  low  temperatures.  One  result  of  this 
refinement  is  a  clearer  picture  of  the  dependence  of  electron  hopping  activation  energy  on 
crystal-field  parameters.  A  generalization  of  the  earlier  exponential  “two-fluid”  model  is  intro¬ 
duced  to  describe  more  accurately  the  high  rate  of  change  in  supercarrier  density  with  tempera¬ 
tures  in  the  regime  where  T  approaches  Tc.  This  new  model  can  provide  additional  insights 
from  the  fitting  of  theory  to  data.  At  the  lowest  temperatures,  supercarrier  densities  may  be 
approximated  by  a  linear  function  of  temperature  if  Tc  values  are  large.  For  lower- Tc  super¬ 
conductors  the  dependence  is  more  typically  exponential.  With  these  improvements,  precise 
fits  to  penetration  depth  versus  temperature  measurements  for  high-purity  YBCO  single  crys¬ 
tals  were  obtained,  as  well  as  more  accurate  agreement  with  critical  (transport)  current  density 
versus  temperature  data  for  thin  films.  Projections  of  critical  current  densities  as  a  function 
of  reduced  temperature  t  =  T!Tc  are  offered  for  the  highest  Tc  materials  predicted  by  this 
theoiy.  Parameter  values  desired  for  optimum  superconductor  properties  are  reexamined. 
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1.  INTRODUCTION 


From  a  study  of  the  origin  of  coherent  charge  transfer  originating  from  large  polaron  conduction 
in  layered  cuprates  at  low  temperatures,  a  theory  of  superconductivity  based  on  covalent  electron  transfer 
(CET)  was  developed  [1,2].  The  fundamental  concept  that  was  used  to  determine  the  transition  tempera¬ 
ture  was  the  statistical  partitioning  of  carriers  between  polaronic  (adiabatic  and  coherent)  and  thermally 
activated  (incoherent)  hopping.  Where  the  electronic  bandwidth  is  much  geater  than  the  polaron  trap 
energy,  i.e.,  b  »  £hop’  it  is  assumed  that  all  carriers  are  large  polarons  that  can  undergo  spontaneous 
delocalization  over  distances  of  several  lattice  parameters,  unless  activated  to  hopping  by  a  thermal 
energy  transfer.  In  the  p-type  cuprates,  charge  transfer  occurs  between  partially  filled  3d  shell  Cu3+  (dz) 
ions  and  the  host  Cu2+  (d  9)  ions;  b  is  determined  principally  by  the  covalent  coupling  of  the  Cu  3d£._yl  orbital 
lobes  with  those  of  the  2 p  oxygen  states.  In  cases  where  n-type  charge  transfer  occurs  between  the  filled 
3d  shell  Cu1+  (d10)  ions  and  the  host  Cu2+  (d9)  ions,  the  natural  overlap  of  the  3d  and  4s  bands  can  permit 
a  larger  contribution  from  the  spherically  symmetric  4s  wave  function.  This  overlap  would  manifest  itself 
as  a  reduced  anisotropy  in  the  basal  planes  of  the  layered  cuprates. 

Polaron  trap  energies  arise  from  electrostatic  fields  of  the  ionic  bond  and  the  crystal  field  interac¬ 
tion  between  cation  and  ligand.  As  determined  from  electron  hopping  activation  energies,  Ehop  is  gener¬ 
ally  on  the  order  of  10 -2  eV  where  spin  polarization  energy  is  not  involved,  e.g.,  mixed  valence  Ti,  Bi, 
and  Tl  compounds.  For  oxides  of  the  3d  transition  metal  series,  antiferromagnetism  can  increase  the 
stabilization  of  the  carrier  in  its  trap  by  at  least  two  orders  of  magnitude.  For  reasons  that  are  related  to 
the  magnetic  dilution  from  zero-spin  ( S  =  0)  Cu3+  and  Cu1+  ions  [3],  the  breakdown  of  long-range 
antiferromagnetic  order  reduces  the  activation  energy  Ehop  to  values  below  10  meV.  Since  the  narrowing 
of  b  with  increasing  temperature  is  not  an  issue  as  long  as  the  large  polaron  condition  b  »  Ehop  is 
satisfied,  the  only  temperature  dependence  of  the  polaron  supercarrier  density  ns  arises  from  a  Boltzmann- 
type  partition  fraction  of  hopping  carriers,  which  become  dominant  where  Ehop!kT  «  1.  Except  for  the 
spatial  ordering  of  the  polaron  supercarriers,  which  is  necessary  to  satisfy  the  London  Vns  =  0  condition, 
this  temperature  dependence  is  analogous  to  that  of  the  collective  carrier  distribution  of  a  degenerate 
semiconductor. 

In  this  report  an  additional  temperature-dependent  factor  is  introduced  to  account  for  the  scattering 
of  polaron  carriers  by  lattice  vibrations.  This  scattering  effect  does  not  reduce  the  density  of  carriers  but 
disrupts  the  formation  of  coherent  chains  of  polaronic  carriers  that  are  the  basis  of  the  superconducting 
state.  Moreover,  the  two-fluid  function  developed  previously  [2,4]  is  placed  in  a  more  general  format  that 
makes  it  sensitive  to  specimen  purity  and  homogeneity. 
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2.  THEORETICAL  DETERMINATION  OF  Tc 

For  a  system  of  N  carriers  per  unit  volume,  the  fraction  of  supercarriers  ns  is  determined  by  a  series 
of  probabilities  [2]:  (a)  that  of  polaronic  conduction,  i.e.,  not  thermally  activated  (hopping)  Ppol  = 
[1  -exp(-£hop/k7)],  (b)  that  of  polaron  isolation  Asol  ”  1  -  A  fix,  where  0  <  p  <  1  is  a  polaron  dispersal 
parameter,  and  (c)  the  quantum  mechanical  transfer  probability  Pttans  =  ??  derived  from  the  overlap  of 
3d^  and  2 p  wave  functions  (p-type  case).  The  resulting  relation  became: 

ns  _  D  p  p 

N  ~~  MxansMsoKpol 

=  77(1  -  4j8r)  1  -  exp^-§^^) 

The  original  expression  for  t]  was  derived  from  a  molecular-orbital  model  that  included  a 
Cu2+-02_-Cu2+  superexchange  calculation,  from  which: 

rj  =  1  —  Q  ,  (2) 

where  W  =  A E/b,  with  b  representing  the  exchange  integral  energy  derived  from  the  molecular-orbital 
theory  of  covalent  bonding  and  A E  =  AE’q'F  describing  the  electrostatic  energy  of  the  polaron  trap.  AE() 
=  (e2/Ka),  where  K  is  the  effective  dielectric  constant,  a  is  the  minimum  separation  between  Cu  sites  in 
the  x-y  plane  (lattice  parameter),  and  is  a  function  of  carrier  concentration  x  that  has  been  modeled 
in  Dionne  [1,2]  to  reflect  four  and  two  interacting  polarons,  respectively.  In  its  simplest  approximation 
*F  =  1  -  2x. 

The  parameter  b  was  originally  chosen  as  0.43  eV  based  on  a  calculation  derived  from  Anderson’s 
theory  of  superexchange  in  transition-metal  oxides  [1,2].  For  the  following  analysis,  it  has  been  found 
appropriate  to  consider  this  quantity  as  a  slowly  decreasing  function  of  polaron  concentration  x.  For  the 
layered  cuprates,  the  following  relation  is  used  to  obtain  Tc  values  that  fit  experiment: 

b  =  b0(l-ax )  ,  (3) 

where  standard  values  of  b{)  =  0.36  eV  and  a  =  0.5  eV  based  on  fits  to  data  for  YBa2Cu307_g  (YBCO). 
Since  b  depends  directly  on  the  overlap  of  orbital  wave  function  lobes  [2],  this  dependence  of  b  on  x  is 
at  least  qualitatively  plausible  because  of  the  anticipated  reduction  in  covalent  linkages  due  to  the  absence 
of  eg  orbital  bond  linkages  that  occurs  when  either  Cu3+  or  Cu1+  is  formed. 

Equation  (2)  will  now  be  viewed  as  an  approximation  to  the  more  general  relation: 

77  =  exp(-Q) 

which  will  be  used  in  future  analysis. 


(la) 

(lb) 


(4) 


It  is  now  appropriate  to  recall  the  carrier  threshold  xs  =  xt  (typically  <  0.1)  for  superconductivity 
[2]  as: 

for  a-0-5 

=  0  for  £2  <0.5  .  (5) 


A  revised  expression  of  Equation  (1)  for  the  supercarrier  fraction  then  becomes: 


N  x 


exp(-Q)(l  -  4/k) 


1  -  exp  - 


(6) 


To  solve  for  Tc  as  worked  out  originally  in  Dionne  [1],  let  xs  =  xt  according  to  Equations  (5)  and 
(6).  With  the  transfer  probability  now  expressed  exponentially,  a  new  set  of  theoretical  curves  of  Tc  as 
a  function  of  x  are  plotted  in  Figure  1  for  the  parameters  in  Table  1.  For  p-type  cuprate  superconductors, 
note  that  op  increases  from  a  relatively  low  value  of  4  meV  to  values  of  12  meV  as  the  strength  of  the 
axial  crystal-field  component  increases  with  the  change  in  the  Cu  environment  from  distorted  octahedral 
to  pyramidal  to  planar.  The  relation  between  the  magnitude  of  ^hop  and  the  oxygen  coordination  about 
the  Cu  ions  is  diagramed  in  Figure  2,  which  depicts  the  increase  in  ^hop  as  a  progressive  destabilization 
of  the  upper  dx2  2  orbital  state  with  increasing  axial  (tetragonal)  crystal-field  strength  A.  These  larger 
values  of  A  have  the  effect  of  further  stabilizing  the  low-spin  (5  =  0)  state  of  Cu3+  ions  [1-3].  Where 
the  transfer  orbital  contains  a  significant  45  wave  function  component,  as  in  the  case  of  n-type  cuprates 
or  metals,  the  influence  of  this  crystal-field  splitting  would  be  correspondingly  reduced.  Consequently, 
Ehgp  and  Tc  values  would  be  smaller.  The  improvement  in  polaron  dispersal  also  serves  to  enhance  Tc 
by  allowing  larger  concentrations  of  polarons  to  participate  in  supercurrents. 
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Figure  1.  CET  theory  projection  of  critical  temperatures  for  cuprate  superconductors  as  functions  of  polaron 
charge  density  and  the  degree  of  polaron  charge  dispersal.  Curve  A  is  the  (LaSr)2Cu04  system  with  distorted 
octahedral  coordinations ,  curve  B  is  YBa2Cuf07  s  with  pyramidal  coordinations,  curve  C  is  the  Tl/Bi/Hg  cuprates 
with  planar  coordinations ,  curve  D  is  curve  C  with  enhanced  b,  and  curve  E  is  a  proposed  theoretical  limit  with 
b  =  0. 
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TABLE  1 

Parameters  for  Calculating  Tc  in  Figure  1 


Curve 

^hop 

(meV) 

P 

K* 

a 

(A) 

bo 

(eV) 

A 

4 

0.5 

16 

4 

0.37 

B 

10 

0.45 

16 

4 

0.36 

C 

12 

0.4 

16 

4 

0.36 

D 

12 

0.4 

16 

4 

0.37* 

E 

12 

0 

16 

4 

0.36 

‘Effective  values  of  K  determined  from  the  stanc 

lard  relation  (VK0  tical 

-MK0)~\ 

For  these  compounds,  K0  = 

25  and  Koptic al  «  10  (estimated). 

"•increased  b0  represents  the  effect  of  a  larger  orbital  overlap  due  to 
increased  pressure  used  in  sample  preparation. 

145391-79 


Figure  2.  Growth  of  the  Cu3+  d8  low-spin  state  (S  =  0)  stabilization  by  eg  doublet  splitting  as  the  tetragonal  crystal 
field  component  progresses  through  distorted  octahedron  (D4h)  pyramid  (C4v)  — >  plane  (C2h)- 


Projections  of  Tc  for  the  hypothetical  case  of  £hop  “  12  meV  and  b  =  0  suggest  that  Tc  may  reach 
values  in  excess  of  300  K.  The  transfer  integral  b,  which  may  be  artificially  increased  by  applying  high 
pressure  to  increase  the  overlap  of  the  covalent  bonding  orbitals,  can  enhance  Tc  of  Tl/Bi-class  com¬ 
pounds  to  levels  above  150  K  by  raising  b0  from  0.36  to  0.37  eV.  Critical  temperatures  in  this  regime 
have  been  reported  for  HgBa2Cas_1CuB02f|+2+g  samples  prepared  under  extreme  pressure  [5]. 
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3.  POLARON  SCATTERING  BY  LATTICE  VIBRATIONS 


Electronic  conduction  in  metal  oxides  with  mixed- valence  transition-metal  ions  is  typically  polaronic 
at  low  temperatures,  where  the  carriers  are  “dressed”  by  local  lattice  adjustments  to  accommodate  the 
change  in  ionic  size  and  electrical  charge  at  the  isolated  carrier  site.  According  to  the  ideas  set  forth  in 
the  Holstein  treatise  [6],  polarons  exist  in  an  ionic  compound  as  a  result  of  the  covalent  bonding 
component  that  broadens  the  electronic  energy  level  of  the  outer  electron  into  a  band  of  width  b.  In 
general,  these  bandwidths  are  smaller  than  the  carrier  trap  or  activation  energy,  i.e.,  b  <  £hop,  so  that  the 
covalent  transfer  or  tunneling  probability  is  proportional  to  the  b/Eho]  ratio.  As  temperature  increases,  the 
phonons  serve  to  stabilize  the  carrier  in  its  trap  by  increasing  its  lifetime,  thereby  reducing  the  influence 
of  the  covalent  interaction.  As  a  result,  the  phonons  tend  to  enhance  carrier  localization  by  having  the 
effect  of  narrowing  the  polaron  band  (which  is  the  electronic  band  narrowed  by  phonons)  and  restricting 
the  polaron  range  to  a  single  lattice  parameter,  i.e.,  creating  a  small  polaron.  Eventually,  the  band¬ 
narrowing  will  remove  the  covalent  influence  and  polaron  activity  altogether,  allowing  incoherent  hop¬ 
ping  to  take  over  as  the  principal  conduction  mechanism  of  these  basically  ionic  compounds. 

For  high-7^  superconductors,  however,  neither  electron  hopping  nor  small  polarons  will  enable  the 
required  energy-free  charge  transport  because  of  the  dependence  on  phonon  interactions.  Large  polarons 
are  necessary  for  conduction  with  a  coherence  length  of  several  lattice  parameters,  and  this  situation 
requires  that  b  »  7?hop.  In  this  case,  isolated  carriers  may  move  spontaneously  over  several  lattice 
parameters,  limited  only  by  the  electrostatic  attraction  between  carrier  and  fixed-site  charge  at  the  impu¬ 
rity  site  or  vacancy  that  created  the  polaron,  i.e.,  the  other  half  of  a  flexible  dipole.  Typically,  b  is  of 
the  same  magnitude  as  £hop  (~10-1  eV)  and  strictly  a  function  of  the  covalent  parts  of  the  chemical 
bonding  of  the  particular  compound.  As  a  consequence,  large  polarons  can  be  created  only  by  reducing 
£hop.  Two  situations  where  large  polarons  are  possible  through  lowering  ■^hop  to  a  range  below  10  2  eV 
by  eliminating  the  magnetic  component  of  the  trap  energy  are: 

1 .  The  spins  of  both  transfer  ions  are  aligned  ferromagnetically  with  each  other  and  with 
those  of  the  surrounding  lattice.  Under  these  conditions,  the  exchange  stabilization  is  the 
result  of  electron  delocalization  through  superexchange  by  covalent  interaction  with  oxy¬ 
gen.  Consequently,  the  stabilization  actually  drives  rather  than  impedes  the  charge 
transfer.  This  process  occurs  with  Mn3+(d4)  <-»  Mn4+(<73)  +  e~  ferromagnetic  couplings 
in  (La3+1__y42+JC)[Mn3+1_JcMn4+JC]03  compounds  that  feature  large  magnetoresistance 
effects  at  the  Curie  temperature.  For  x  <  0.5,  the  ferromagnetism  is  reinforced  because  the 
Mn3+-02”-Mn3+  superexchange  is  also  ferromagnetic  through  the  influence  of  a  quasi¬ 
static  Jahn-Teller  effect.1 


1 .  The  compexities  of  the  various  superexchange  interactions  involving  the  different  valence  states  of  the  Mn  ions 
were  examined  in  J.B.  Goodenough,  A.  Wold,  R.J.  Amott,  and  N.  Menyuk,  “Relationship  between  crystal  symmetry 
and  magnetic  properties  of  ionic  compounds  containing  Mn3+,”  Phys.  Rev.  124,  373  (1961). 
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2.  One  of  the  transfer  ions  has  zero-spin  and  the  lattice  is  magnetically  disordered.  In  this 
case,  there  is  no  internal  exchange  stabilization  gained  from  satisfying  the  Hund’s  rule  spin 
polarization  requirement,  and  neither  ion  has  an  exchange  commitment  to  its  lattice  envi¬ 
ronment.  Therefore,  there  is  no  static  exchange  stabilization  of  any  kind.  This  condition 
is  precisely  that  of  the  Cu2+(<i9)  <H>  Cu3+(<i8)  zero-spin  +  e~  transfer  of  the  superconducting 
layered  cuprates  discussed  in  Dionne  [1,2]. 

Both  situations  are  readily  identifiable  by  the  metallic  nature  of  the  resistivity  at  low  temperatures. 

For  a  material  with  x  carriers  per  formula  unit  of  dimension  a,  separation  between  polaron  sources 
Lp  =  a/x.  If  this  distance  is  treated  as  a  mean  free  path  for  coherent  transfer  of  carriers  between  large 
polaron  cells,  scattering  by  lattice  vibrations  may  decrease  the  density  of  available  supercarriers  accord¬ 
ing  to: 


(  L  \ 

Aw,Sv  ~  ~nSV  ^  'j 


(7) 


Since  1/(AV)  =  (1  /XD)(T/0D),  where  XD  and  ©D  are  the  respective  Debye  wavelength  and  tempera¬ 
ture,  and  0D~  hVj/k  with  the  Debye  frequency  expressed  in  terms  of  the  acoustic  phonon  velocity  as 
VA  =  va  /  XD,  Equation  (7)  may  be  reduced  to: 

= -»4f)(£)te)  ’  (8) 

where  va  ~  a{C/M)m  ~  ( Y/ps)m .  The  variables  involved  in  this  approximation  are  the  lattice  elastic  force 
constant  C,  the  effective  atomic  mass  M,  the  Young’s  modulus  Y,  and  the  mass  density  ps. 

After  solving  the  simple  differential  equation,  a  lattice  scattering  probability  may  be  defined: 


where  L,  =  (k/h)(a/va).  For  the  superconducting  cuprates,  a  ~  4  x  10  8  cm,  va  ~  6.5  x  105  cm/s,  and 
£  =  1.3  x  10-3  K”1.  In  terms  of  the  reduced  temperature: 


Platt  -  exp 


(10) 


The  scattering  probability  would  be  significant  where  Tc  is  large  but  less  important  with  lower- Tc  mate¬ 
rials,  e.g.,  metal  and  intermetallic  compounds. 
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4.  GENERALIZED  TWO-FLUID  FUNCTION 


In  the  foregoing  analysis,  an  expression  for  the  density  of  available  supercarriers  ns  was  calculated 
as  a  function  of  the  total  carrier  density  N.  However,  the  variable  ns  does  not  represent  the  actual  carriers 
participating  in  the  supercurrents  because  only  an  effective  fraction  (ns/N)  of  these  polarons  will  be 
capable  of  such  ordered  linkages  at  any  time  for  at  least  two  reasons.  First,  scattering  by  lattice  vibrations 
will  cause  deflections  or  misdirections  of  the  carriers,  particularly  as  the  temperature  increases  and  the 
average  wavelength  becomes  shorter;  second,  inhomogeneities  arising  from  chemical  impurities  or 
microstructural  inhomogeneities  will  prevent  a  general  onset  of  complete  supercurrent  near  7\  Only  at 
T-  0  may  one  reasonably  assume  that  nse  =  ns.  A  general  function  capable  of  representing  the  growth 
of  supercurrent  near  Tc  may  be  derived  as  follows: 

nse(T)  =  ns(T)-nt(T)  ,  (11) 


where  the  threshold  supercarrier  density  nt(T)  =  nt(Tc)[  1  -  nse(T)/ns( 0)].  Upon  substituting  this  relation 
into  Equation  (11),  the  result  is: 


p  _  nse(T)  _  ns(T) 
rl  “  «,(0)  ns(0) 


n,(Tc ) 
ns(0) 


[1-Fj] 


(12) 


or 


p  _  ns(T)-n,(Tc) 
rl  "  ns(0)-n,(Tc) 

Based  on  the  logic  of  Equation  (12)  the  general  relation  may  be  written  as: 


or 


p  _  ns(F)  nt(Tc)  i-i  p  -on 

Fm-W0)— 


1  yyi  —  ..  /n\  „  /n\  J 


(13) 


(14a) 


(14b) 


m~ns(  0)  ns(  0)' 

where  t  =  T/Tc.  If  W  =  Ehop/kTc,  the  simple  two-fluid  relation  of  Equation  (12)  may  be  expressed  as: 


Fx  =l-exp(w--^) 


(15) 


which  is  the  original  relation  from  the  CET  model  [2,4],  The  general  version  of  Equation  (14)  becomes: 


Fm  =  1  -  exp(-f)  -  [1  -  exp(-W)][exp(  W  -  f] 


W-,m 


(16) 


To  illustrate  the  influence  of  the  parameter  m,  Equation  (16)  is  plotted  in  Figure  3  as  a  function 
of  t  with  W  =  1.25  for  values  of  m  ranging  from  1  to  100.  In  the  limiting  case  m  — »  °°,  the  onset  of  the 
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superconducting  state  would  be  sudden  and  complete  as  t  — >  1 .  This  situation  would  be  ideal  and  require 
perfect  structural  homogeneity  and  optimum  dispersal  of  polaron  sources.  Therefore,  m  can  be  treated  as 

a  material  quality  parameter. 
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Figure  3.  Two-fluid  function  Fm  for  W  -  1.25  and  m  =  1,  2,  4,  8,  and  100. 

When  the  lattice  scattering  probability  factor  Platt  of  Equation  (10)  is  included,  the  carrier  variables 
are  redefined  as: 

ns(t)  =  N  exp(-^j  exp(-Q)(l-4/k)[l-exp(-J f)] 

ns(0)  =  N  exp(-Q)(l  -  4 fix)  ,  (17) 

ns(l)-N  exp^--^^  exp(-Q)(l  -  4/k)[l  -  exp(-W)] 
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For  the  YBCO  case,  x  =  0.25,  Tc  =  93  K,  W=  1.25,  and  Equation  (10)  reduces  to  Platt  =  exp(-t/2).  The 
general  two-fluid  function  of  Equation  (14b)  is  plotted  in  Figure  4  for  the  same  range  of  m  values  with 
the  data  of  Bonn  et  al.  [8].  In  Figure  5,  the  inverse  ratio  [nse(0)/nse(t)]m  =  Fm~m  (which  represents  the 
temperature  variation  of  the  London  penetration  depth)  is  plotted  for  the  same  set  of  parameters,  except 
that  Tc  is  reduced  slightly  to  90  K  and  compared  with  the  data  selected  from  Mao  et  al.  [9].  For  m  =  8  in 
both  cases,  excellent  agreement  was  obtained  with  the  data  that  were  derived  from  microwave  measure¬ 
ments  of  penetration  depth  on  nearly  perfect  single  crystal  (s.c.)  specimens.  Figure  6  presents  an  inter¬ 
esting  comparison  between  theory  and  experiment  for  a  YBCO  textured  film.  From  the  data  of  Porch 
*  et  al.  [10],  agreement  is  established  for  m  =  4  at  larger  t  values,  but  a  growing  departure  is  apparent  at 

low  t  values.  The  failure  to  produce  a  closer  fit  at  low  t  could  be  accounted  for  by  a  depletion  in  the 
long-wavelength  phonon  populations  due  to  the  presence  of  grain  boundaries  along  longitudinal  paths, 
'  as  well  as  the  submicrometer-scale  thickness  of  the  film.  The  absence  of  linearity  in  the  slope  near  t  =  0  could 

result,  for  example,  from  an  inverse  quadratic  relation  between  (A.v)  and  T,  i.e.,  H(kv)  ~  (lfkD)(T/®D)2  where 
T  «  ®jy  thereby  introducing  a  T2  dependence  to  the  exponent  of  Equation  (9).  The  parameters  used  in 
these  computed  temperature-dependent  properties  and  those  to  follow  are  summarized  in  Table  2. 
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Figure  4.  Two-fluid  function  Fm  with  phonon  scattering  included  for  W  =  1.25  and  m  =  1,  2,  4,  8 ,  and  100.  In  this 
case,  the  YBCO  parameters  of  x  =  0.25  and  Tc  =  93  K  are  used.  Single-crystal  data  are  from  Bonn  et  al  [8]. 
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2.0 


Figure  5.  Penetration  depth  ratio  [nse(0)/nse(t)]l/2  (=  Fm~I/2)  as  a  function  oft  compared  with  s.c.  data  of  Mao 
et  al.  [9]  for  YBCO  with  W  =  1.29  and  Tc  =  90  K.  Other  parameter  values  used  are  the  same  as  those  for  Figure  4. 


W 
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TABLE  2 

Parameters  for  Calculations  in  Figures  2  through  8 


Sample 

Type 

m 

y* 

(x  10'3) 

^hop 

(meV) 

m 

W 

Ref 

YBCO 

crystal 

8 

5.2 

10 

93 

1.25 

8 

YBCO 

crystal 

8 

5.2 

10 

90 

1.29 

9 

YBCO 

text,  film 

4 

5.2 

10 

90 

1.29 

10 

YBCO 

text,  film 

3 

5.2 

10 

87 

1.33 

11,12 

NbN 

gran,  film 

2 

5.2* 

2.7* 

11.9 

1.76* 

13 

TBCCO 

text,  film 

3 

5.2 

12 

105 

1.32 

14 

(TI2Ba2CaCu208) 

‘Undoubtedly,  x  is  larger  for  this  compound,  but  £  is 
phonon  velocity. 

probably  also  larger  because  of  a  smaller 

tValue  consistent  with  W  =  1 .76  and  Tc 

=  1 1 .9  K. 

tValue  based  on  typical  Bardeen-Cooper-Schrieffer  (BCS)  gap  energies. 

Another  important  property  of  superconductors  that  may  be  interpreted  with  this  model  is  the 
critical  current  density  ic  ,  usually  measured  as  a  transport  current  and  divided  by  the  cross-sectional  area 
of  the  specimen  or  film.  In  the  most  recent  study  of  CET  [2],  i£  was  related  to  the  effective  supercarrier 
density  of  a  thin  film  according  to  n se312.  If  this  model  is  combined  with  the  preceding  results  for  nse(t), 
the  following  is  obtained: 


UP  _  F  3/2 
ic(0)  rm 


(18) 


In  Figure  7,  Equation  (18)  is  plotted  for  thin  film  YBa^UjOy,  and  data  used  in  earlier  work  [1 1,12] 
are  replotted  to  show  a  substantially  improved  fit  with  theory,  in  this  case  with  a  lower  value  of  m  =  3 
that  appears  to  produce  a  better  fit  to  textured  film  specimen  data.  For  an  additional  test  of  the  theory, 
the  case  of  the  BCS  superconductor  NbN  (Tc  =  1 1.9  K)  is  plotted  for  W  =  1.76,  and  m  =  2  (granular  film) 
together  with  experimental  points  of  Clem  et  al.  [13].  Another  set  of  data  [14]  from  commercial  TBCCO 
with  Tc  =  105  K  plotted  in  logarithmic  form  in  Figure  8  compares  excellently  with  theory. 


v 
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Figure  7.  Reduced  critical  current  density  ic(t)/ic(0)  as  a  function  of  t  for  YBCO  thin  films,  with  W  =  1.33  and  Tc  = 
87  K.  Data  are  from  Westerheim  et  al  [11]  and  Inam  et  al  [12].  The  second  calculated  curve  is  for  NbN  with 
W  =  1.76 ,  the  BCS  value ,  compared  with  data  of  Clem  et  al.  [13].  For  the  oriented  films  of  YBCO,  m  =  3;  for  the 
granular  film  of  NbN,  m  is  reduced  to  2. 


Another  exercise  of  interest  is  to  apply  this  new  model  to  a  point  on  the  (proposed)  optimum  case 
curve  E  of  Figure  1  to  estimate  the  effects  of  the  likely  trade-offs  between  Tc  and  ic  as  the  critical 
temperature  increases.  To  this  end  Equation  (18)  is  plotted  in  Figure  9  for  parameter  values  of  j3  =  0,  x 
=  0.25,  and  m  =  8.  Here  W  is  reduced  to  0.46  because  Tc  is  raised  to  300  K,  and  ^hop  remains  at  12  meV. 
The  reduced  ic  curve  for  s.c.  YBCO  with  Tc  =  93  K  and  W  =  1.25  is  plotted  to  illustrate  the  advantage 
in  current  capability  at  T  =  77  K,  where  ic  is  more  than  three  times  greater  for  the  300-K  superconductor. 
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Figure  9.  Comparison  of  ic(t)/ic(0)  versus  t  for  standard  YBCO  with  a  postulated  300-K  cuprate  superconduc 
tor.  In  both  cases ,  m  =  8y  the  value  for  high  quality  single  crystals. 


5.  DISCUSSION 


Refinements  to  the  CET  model  for  determining  Tc  indicate  that  hopping  electron  activation  energies 
should  be  more  than  double  (>10  meV)  those  reported  in  the  original  work  [1,2]  and  are  directly  depen¬ 
dent  on  the  strength  of  the  axial  crystal-field  component  around  the  Cu  ions.  As  concluded  in  the  initial 
work  on  this  theoretical  approach  [1],  critical  temperatures  could  exceed  300  K  through  the  attainment 
of  ideal  dispersal  Q 3  — >  0)  and  larger  activation  energies.  Reports  of  superconductivity  behavior  at 
temperatures  near  250  K  [15]  in  cuprates  prepared  by  a  sequential  epitaxial  layer  growth  technique 
support  the  earlier  suggestion  that  large  ^hop  values  obtained  from  planar  oxygen  coordinations  and 
nearly  perfect  cation  or  vacancy  dispersals  by  some  type  of  individual  layer  deposition  method  could  lead 
to  optimized  cuprate  superconductor  structures. 

The  influence  of  phonon  interactions  on  polaron  charge  transfer  is  characterized  as  an  exponential 
factor  that  causes  an  additional  decrease  in  the  effective  supercarrier  density  with  increasing  temperature. 
The  simple  two-fluid  model  derived  from  this  theory  has  also  been  generalized  to  describe  variations  in 
the  growth  of  supercurrent  near  T  .  One  notable  result  of  this  analysis  is  that  the  variation  of  supercarrier 
density  with  temperature  becomes  linear  at  temperatures  near  zero  if  Tc  is  large  and  the  material  is  in  the 
form  of  a  high  quality  single  crystal.  For  textured  films  there  is  evidence  that  the  linear  relation  does  not 
hold  at  low  temperatures,  possibly  because  of  a  nonlinear  relation  between  the  mean  phonon  wavelength 
and  temperature  due  to  the  presence  of  grain  boundaries  and  submicrometer  film  thicknesses.  For  lower- 
Tc  materials,  metals  for  example,  the  temperature  dependence  returns  to  the  exponential  behavior  pre¬ 
dicted  by  the  earlier  model  without  lattice  scattering  effects  or  that  typical  of  the  BCS  theory.  In  these 
cases,  the  two-fluid  universal  curves  would  resemble  those  of  Figure  3  rather  than  Figure  4  at  lower 
values  of  t. 

As  a  consequence,  higher- Tc  materials  would  have  supercurrent  capacities  that  are  more  vulnerable 
to  lattice  vibrations  near  their  highest  operating  temperatures.  From  the  calculated  curves  plotted  in 
Figure  9,  the  relative  supercarrier  density  at  90%  of  Tc  for  a  material  with  Tc  =  300  K  is  only  a  small 
fraction  of  that  for  YBCO  with  Tc  =  93  K.  Nonetheless,  there  is  a  substantial  advantage  to  higher  critical 
temperatures  for  operation  in  the  lower  temperature  regimes.  A  300-K  superconductor  engineered  from 
these  systems,  for  example,  theoretically  would  have  almost  a  fourfold  advantage  in  ic  over  YBCO  if 
operated  at  77  K.  If  a  system  could  be  found  with  a  greater  ^hop’  the  W  parameter,  hence  the  ic(t)  curve, 
would  be  correspondingly  higher. 

i  With  these  revisions  to  the  CET  theory,  the  parameter  trends  for  the  design  of  a  hypothetical 

superconductor  with  optimum  properties  are  reexamined.  To  achieve  the  highest  T„  b  should  be  as  large 
as  possible  and  £hop  as  large  as  tolerable  (bearing  in  mind  the  implications  of  magnetic  exchange  trap 

4  contributions  to  £hop  and  the  spin-alignment  frustration  requirement  discussed  in  Dionne  [2]).  These 

conditions  will  permit  a  corresponding  increase  in  A E  to  preserve  a  reasonably  small  Q  (=  A E/b)  and 
maintain  a  low  threshold  polaron  density  xf  (In  low-7],  metals,  for  which  b  is  small  because  of  the  lower 
ionization  potentials  of  the  outer  electrons,  xt  can  be  large  [2]).  If  A E  is  raised  through  a  decrease  in 
temperature,  there  is  an  additional  benefit  of  an  increase  in  the  condensation  energy  AG  as  determined 
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by  the  ferroelectric  alignment  of  polarons  [2].  With  a  larger  AG,  higher  critical  magnetic  field  and  critical 
current  densities  would  result.  Finally,  to  reduce  the  effect  of  lattice  vibration  scattering  of  coherent 
polarons,  an  increased  acoustic  phonon  velocity  is  desired,  requiring  high  stiffness  and  low  mass  density. 

The  layered  cuprates  are  only  partially  covalent,  and  their  dielectric  constants,  although  lower  than 
those  of  metals  and  intermetallics,  are  still  relatively  high,  particularly  in  the  larger  cation  compounds  in 
which  superconductivity  has  been  found.  There  are,  however,  other  possible  vehicles  for  large  polaron 
superconductivity.  Recent  work  in  polymers  [16]  may  provide  advances  in  this  field,  particularly  in  view 
of  their  high  covalency,  low  dielectric  constants,  and  lower  mass  densities. 
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densities  may  be  approximated  by  a  linear  function  of  temperature  if  Tc  values  are  large.  For  lower- Tc  supercon¬ 
ductors  the  dependence  is  more  typically  exponential.  With  these  improvements,  precise  fits  to  penetration  depth 
versus  temperature  measurements  for  high  purity  YBCO  single  crystals  were  obtained,  as  well  as  more  accurate 
agreement  with  critical  (transport)  current  density  versus  temperature  data  for  thin  films.  Projections  of  critical 
current  densities  as  a  function  of  reduced  temperature  t-T/T  are  offered  for  the  highest  Tc  materials  predicted  by 
this  theory.  Parameter  values  desired  for  optimum  superconductor  properties  are  reexamined. 
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